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Using the idea of Itzykson-Zuber integral, unitary-matrix integration of 2D 
Yang-Mills action is presented. The uniqueness of the solution of heat equation 
enables us to integrate out the unitary-matrix parts of hermite matrices and to 
obtain the expression of integration over vectors, also in this case. 
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Introduction 

In 1980, C. Itzykson and J.-B. Zuber presented the method of integration over angular 
variables on a simplest model coupling two hermite matrices pQ. Here, the angular variable 
means the unitary matrix which diagonalizes physical quantitiy, i.e. NxN hermite matrix. 
Their work had a great impact and enabled developments especially on the theories of 
random matrix and 2D gravity ( etc.). While, until then, there was Creutz's unitary- 
matrix integration |3] well known in lattice gauge theory which integrate unitary matrices 
order by order in high temperature expansion, Itzykson-Zuber integral is epoch-making in 
the sense that they obtained the expression only using the eigenvalues of hermite matrices, 
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after exactly integrating out angular variables, to all orders of exponential. In practice, 
since it is too complicated to integrate in usual ways, they evaluate the integral with the 
partial differential equation known as heat equation whose uniqueness of its solution is 
completely understood. They showed that the expression after integration is equivalent 
to that before integration using the uniqueness of the solution. After the unitary-matrix 
integration, the number of variables is reduced to N, i.e. vector, from NxN matrix, and 
therefore, the analyses such as N—>oo planar limit becomes very easy. 

In this paper, using their idea, we present a similar expression after integrating out 
unitary matrices on 2-dimensional Yang-Mills action with intent to prepare to apply this, 
in future, to 10-dimensional Yang-Mills theory, i.e. IIB Matrix Model [HE] which is known 
to be convergent [6J. On our model, if we introduce an external field as NxN hermite 
matrix, the same heat equation as Itzykson-Zuber integral is also found to be satisfied, 
and therefore, we can evaluate our integral on the basis of the uniqueness of its solution. 
However, unlike Itzykson-Zuber integral, since a number of coefficients which cannot be 
determined only from the uniqueness remain, we derive them from the structure of the 
commutator characteristic of Yang-Mills action. As a result, in case N is even, after 
integrating out angular variables, we find that the partition function of 2-dimensional 
Yang-Mills action can be expressed in terms of the remaining eivenvalues. 



Unitary-Matrix Integration of 2D Yang-Mills Action 

The integral we investigate is the Yang-Mills action in 2 "space-time" dimensions 



which exists in the "external field" C. Here, A, B and C are NxN hermitian matrices, 
and iC[A, B] is also hermitian. dA (and dB) is the measure 





N N 



dA = Y[ dAa Y[ d(ReAij)d(ImAij), 



(2) 



i=l i<j 

which is invariant under the adjoint operation of U(N) matrices 



diy^AV) = dA 



V e U{N). 



This measure can be written in the diagonalized form as follows: 
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(3) 
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where A is diagonalized by unitary matrix U such as A = Wall, a = diag{a,\, • • • , ajv), 
and dU is the normalized U(N) Haar measure, fdU = 1. A(a) is the Vandermonde 
determinant 
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As we can check easily, the Itzykson-Zuber integral pQ, 

K(C; g) = J dXexp{-g trX 2 + trCx} 



d N x 



A{x) 



exp < — g\x\ + c ■ x 



(5) 



obeys the following partial differenial equation, 
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dC 2 
( d 2 \ 
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-|/f(C; 9 ), 



(6) 



. . , dCij dCji 



Essentially, this equation is the heat equation, and has an unique solution (5) under the 
boundary condition K(C; g— >oo) =0. 

Similarly, the integral (1) satisfies the same equation 



d_ 

dg 



I(C; 



(7) 



and must have an unique solution under the boundary condition I(C; g^oo) = 0. After 
integrating over U(N) matrices, we expect that (1) becomes the integral which only 
contains the integration over two eigenvalues of hermite matrices A and B. Then, we 
suppose the integral (1) to become like the following form: 



i(c;g) = J d N ad N bexp{-g ^ Ri jk iaibja k bi + ^ r ijfc qaA}- (8) 

i,j,k,l=l i,j,k=l 

Here, the coefficients Rijki and are independent of c= (ci, • • • , cat). Since this integral 
must obey the partial differential equation 
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a d 
J2d^ j i ^9) = —^^9), (9) 



Rijki is determined from T^: 



N 



Rijkl ^ ^ ^mij^mkl ■ (^) 



m=l 



Therefore, it is sufficient we obtain r ijfe , and, from the equation 



f)2 N rj2 

3=1 i 

C = UlcUc, 

where V^ c only depends on Uc, the integral (1) becomes 



N N 



d N ad N b exp j - g R ijkl aibja k bi + r^a,-^}, (12) 

i,j,k,l=l i,j,k=l 

owing to the uniqueness of the solution of (7). Here, f(a, b) is the undetermined function 
of a and b. Furthermore, though c is exactly the eigenvalue of hermite matrix C, we 
cannot certain that both a and b are those of A and B, yet. However, with g = in (12), 
the integral 

dl/dV exp 

= J dUdVdW exp{itrW^CW[tfAU,V^BV}} 
dVdW exp |? irA^BV, ^C^]} 
dUdW exp trS^CH 7 ", [/Ut/]} 

N 

exp | ryfeQaj6 fe | (13) 



A(a)A(6)A(c) 



i,j,k=l 



enables us to derive /(a, 6) = A(a) A(6) and make sure that both a and b are the eigenvalues 
of A and £>, because all A, B and C stand on the same ground. And, also from the facts 
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that trC[A,B] = trA[B,C] = trB[C, A] and trC[A, B] = -trC[B, A], the index of Y ijk 
has cyclic symmetry: = Yj ki = ^kij and anti-symmetry under the interchange of its 
indices: T^k — — Y ik j. Therefore, we obtain 

J(C; g ) = J dAdB exp | g tr[A, B] 2 + i trC[A, B] } 

/A(a)A(b) N N 

d N ad N b — — expj-# ^ Rijkiaibja k bi + r ijfc Q%&fc}. (14) 

i,j,k,l=l i,j,k=l 



The remaining task we must do is to determine the coefficients Y^. Since to have 
the matrices A and B multiplied by the constants k is equivalent to make the eigenvalues 
ka and kb respectively, Y^ are also independent of a and b. Then, we can evaluate 
through the structure peculiar to the commutator [*,*]. We can easily check the integral 

/(a, b;g) = J dU A dU B exp {g tr[U{aU A , U B bU B \ 2 } (15) 

N 

<>x p{ _ 5 R ijkiaibja k b^ (16) 

i,j,k,l=l 



follows from (14), and (15) satisfies the differential equation due to the characteristic 
structure of the commutator: 



N r, 

J2 7 -I(d,b;g) = 0, (17) 
i=i OCli 

therefore, from (16), we obtain the algebraic equation which Rijki must obey: 

N N 

RijkibjCikbi = ^ Y mi jY m kibja k bi = 0. (18) 

i,j,k,l=l i,j,k,l,m=l 

This is the identity about and 6j. By introducing d m = i F m ki a kbi, this reduces to 
the identity 

N 

^ ^mijbjd m = (19) 
i,j,m=l 

about bj and d m , and we obtain conditions for Y^: 



5 



N 

^ijk = 0, for fixed j, k. (20) 

i=i 

When N = 4, these equations become 



1^312 


+ 1^412 


= o, 


Tl23 


+ 1^423 


= o, 


1^213 


+ r 4 i 3 


= o, 


Tl24 


+ 1^324 


= 0, 


1^214 


+ r 3 i4 


= o, 


Tl34 


+ T234 


= 0, 



then, except overall factor, all T^s are determined exactly as 

Tl23 = — Tl24 = Tl34 = — 1^234 = 1- (21) 

When N = 2 and N = 3, = for all i,j, k. When N > 4, from the symmetry of its 
indices, has jyC 3 independent components, but (20) are atC 2 equations at most. So, 
we cannot determine all from (20) in general. However, as is obvious from (14), all 
cii,bi and q are contained in the same weight. That is, except the overall factor which 
can be absorbed into the integration measure, we can set | r ij7 - fe | = 1. Then, for even iV>4, 
we can find the solution of (20): 

r ijk = (-iy +j+ \ for i<j<k. (22) 

When N is odd, the solution such as |r^- fc | = 1 does not exist, because odd numbers of 
±1 cannot add up to zero. In fact, when iV = 5, 



1^312 
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+ r 5 i2 
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+ 1^524 


= 


1^213 


+ r4i3 


+ r 5 i 3 


= 0, 


Tl25 


+ 1^325 
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r2i4 


+ r3u 
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= 0, 


Tl34 


+ 1^234 


+ T534 
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+ r 3 i 5 


+ r4i5 


= o, 


Tl35 


+ 1^235 


+ 1^435 
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Tl23 


+ 1^423 


+ 1^523 


= o, 


Tl45 


+ 1^245 


+ 1^345 


= 



the solution does not exist. This strange feature that we can use the expression (14) 
only for even iV seems to come from the following characteristic of the commutator [*, *]. 
The algebra of hermite matrices has the subalgebra which consists of real symmetric 
matrices = A T = A, and, in this subalgebra, the commutator [A, B] always becomes the 
anti-symmetric matrix [A, B] T = — [A, B] which has even rank. That is, the commutator 
[*,*] essentially contains the structure of even rank which is hidden by larger U(N) 
symmetry before the integration, and then, the integration of U (N) matrices may reveal 
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this structure. However, since we need the integral after taking the limit N — > oo, it is 
sufficient for our purpose to have the result that N is even number. 
Lastly, let us line up the result we obtain in this paper. 




dAdB exp <^ g tr[A, B} 2 + i trC[A, B] 




N 

Rijkl ^ r 'mij-T 
m=l 

p r r 

L ijk L kij 1 



mkl j 



ikj ) 



for i < j < k. 



Discussions 

Using the uniqueness of the solution for heat equation which plays a crucial role 
on Itzykson-Zuber integral, we derive an expression which contains eigenvalues of her- 
mite matrices only as integration variables after integrationg out angular variables on 
2-dimensional Yang-Mills action. Then, characteristic structures of the commutator are 
important and needed to determine coefficients in the expression. And the expression is 
possible only for even N, reflecting the structures. 

Our expression can be applied to various problems. For example, as the model of 
Itzykson-Zuber integral which can be solved exactly in iV — > oo planar limit, we can 
evaluate similarly on 2-dimensional Yang-Mills model. And, as refered in introduction, 
by extending our result to larger dimensions, for example IIB Matrix Model, we may 
be able to treat the phenomenon of the phase-transition of space-time, i.e. inflation. 
However, when we extend to higher dimensions than three, we cannot probably use the 
uniqueness of the solution of heat equation, so must introduce other methods. 

Furthermore, since Vandermonde determinant A(x), which always emerges in case 
of integrating over hermite matrices, is the wave function of ground state of Calogero- 
Sutherland model describing the system of fermions on circle 0, and the orthogonal 
polynomials containing A(x) satisfy particular recursion relations dj, we may be able to 
calculate exactly over remaining eigenvalues. 
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